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I. INTRODUCTION 

In recent years, quantum creep of vortices in high T c -superconductors has attracted considerable interest as, on 
the one hand, this phenomenon is responsible for the dissipation of energy and thus is relevant from a technological 
point of view, while, on the other hand, it represents an interesting example of a macroscopic quantum phenomenon, 
the theoretical study of which is challenging. A particularly well defined and technologically relevant configuration 
are vortices trapped by columnar defects introduced into the sample by heavy ion irradiation. Measurements of the 
critical current < 
vortex dynamic; 



nsity and the magnetic relaxation rate show the strong influence these pinning centers have on the 
. In this paper we consider the problem of the depinning of a pancake vortex governed by Hall 
dynamics from a columnar defect in a layered superconductor in the presence of a small (J -C j c ) transport current. 
The external magnetic field is chosen parallel to the c-axis of the superconductor. In the limit j — > the problem 
is semiclassical and the Euclidean action can be calculated in the whole temperature range. In the limit T = 
the problem discussed above has been considered in RcfEl, however, it appears that the approximations made are 
too rough, leading to an inexact result for the decay rate. The main goal of the present work is to improve on the 
analysis of Ref.El and to extend the calculations to the entire temperature range. We adopt the semiclassical approach 
instead of the lowest Landau level (LLL) approximation used in Ref.cl. For^ general review concerning the decay of 
metastable states see the review of Hanggi, Talkner, and Borkovec, see RefJ3. The problem of quantum and classical 
Hall creep of vortices in waripus geometries and for different driving forces F < F c and F <C F c has been studied by 
various authors, see RefsETt2l. 

The process of quantum tunneling is described by a time-dependent saddle-point solution. Consequently, the 
calculations of the decay rate require the specification of the vortex dynamics. The different contributions to the 
dynamics considered traditionally are of the massive, the dissipative, and the Hall type. The (low frequency) equation 
of motion of a single vortex can be written in the form 

<!>o 

— j An — VUpin = —r]v + mv + qv A n, (1) 

where the dynamical forces are balanced by the Lorentz and pinning forces. In conventional superconductors one can 
neglect the contribution of the mass and Hail terms, whereas in high-T c superconductors the Hall force may become 
relevant. In particular, it is widely believedtJ that at low temperatures the superclean limit can be reached where the 
Hall term is large. In this case Wqt jSj 1, where lo q is the leveLspacing inside the vortex core and r is the quasiparticle 



relaxation time. Indeed, recent Hall angle measurements! 1 2 H 13 I demonstrate that the limit a > rj can be realized. On 
the other hand, the large parameter ujqt also gives rise to a large vortex mass. Microscopic, calculations show that 
the vortex mass is enhanced by a factor (ep/A) 2 ~ 100 in comparison to the dirty limitliiTO. Still, it can be shown 
that for frequencies lu < luq the Hall force wins over the inertial one, whereas at high temperatures w > u>q the vortex 
equation of motion cannot be cast into the simple form (jj]) (an accurate description produces dispersive transport 
coefficients f?(w), and a(u)). In this situation a good starting point is to ignore the contribution of the vortex mass and 
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solve the remaining Hall tunneling problem. The advantage of this treatment lies in the fact that in the limit j <C j c 
the problem allows for an analytical solution in the whole temperature range with a universal answer: the Euclidean 
action depends only on the depth of the pinning potential, the detailed shape of the potential being irrelevant. In a 
second step we establish the consistency of this approximation in the physically relevant regime of parameters. The 
outline of the paper is as follows: In section [n] we discuss the model and the qualitative picture of the tunneling 
process. In section III we calculate the Euclidean action in the whole temperature range and show that the problem 
always exhibits a second-order like transition from quantum to classical behaviour. Furthermore, we provide estimates 
for the preexponential factors in the various regimes. Finally, in section [V we discuss the conditions of applicability 
of our results. 



II. MODEL AND QUALITATIVE PICTURE 



Consider a pancake vortex with a Hall-force dominated dynamics trapped in a 2D potential well Uq 
and subject to an external force F, i.e., the effective potential U (x,y) takes the form 



U(x,y)=U (y X 2 + y^ 



Fx. 



(2) 



The function Uo (r) is supposed to be monotonously increasing, Uq (0) = and Uo (oo) = Uo- At distances r much 
larger than the characteristic radius a of the pinning potential but still smaller than the magnetic-field penetration 
length X a b in the a6-plane, Uo (r) behaves as Uo — B/r 2 with B a constant of order Uoa 2 (see Ref.liH). 




Fig.l Inverted potential U after the transformation t — + —ir, y — > iy. Note that the inversion affects only the unstable direction 
along the a;-axis. The equipotential lines define the quasiclassical trajectories. The thick solid line corresponds to the zero temperature 
instanton, whereas the thick dotted line is a finite temperature bounce trajectory. 



The present problem is equivalent to that of the motion of a charged particle in a strong magnetic field. In this 
case, the appropriate Lagrangian takes the form 



L = axy — U (x, y). 



(3) 
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The precise conditions allowing us to neglect the mass of the particle will be discussed in section IV. The particle 
whose motion is described by the Lagrangian (j^) moves along the equipotential lines of the potential landscape U(x,y). 
The Euclidean action of the particle as obtained through the substitution S — J Ldt — > —iS and t — * —it can be 
written in the form 

r+h/2T 

Se U c1 = / [-iaxy + U (x, y)] dr. (4) 

J-H/2T 

Obviously, the imaginary unit appears in the Euclidean action, i.e., the saddle-point solution-|is-,in general complex. 
However, in the case studied here it is possible to reduce the complex problem to a real onecHlil via performing the 
additional transformation y — > iy. If the potential U(x,y) satisfies the condition lm{U(x,iy)} = 0, we obtain a 
real-time problem. In addition, if after the y — > iy-transformation the Lagrangian (Q) exhibits a saddle-point solution, 
we can find it and calculate the decay rate. One can easily see from Eq. (|2j) that the condition Im {U(x, iy)} = 
is indeed satisfied for our potentialjLia i.e., we can study the effective problem of the tunneling of a particle whose 
dynamics is described by the Euclidean action 

+n/2T 



S Eucl [x(T),y(T)}= [aiy + U ( y/x 2 - y 2 )-Fx] dr. (5) 

J-H/2T 

In Fig. 1 we show the inverted potential after the transformations t — > — vr and y — > iy together with the equipotential 
lines. Obviously, this construction produces a new potential shape exhibiting a bounce solution (saddle-point). Note 
that the inversion affects the potential only in the unstable direction along the force. 

As the potential energy U {x, iy) is preserved during motion, the imaginary time trajectories satisfy the equation 

Uq f x 2 — y 2 ^ — Fx = Const. In Fig. 2 we plot these trajectories for the problem at hand. At zero temperature and 



in the limit F — » wc obtain Uq y\/ x 2 — y 2 j — Fx = 0, i.e., SeucI — a § ydx = aA, where A is the area encircled 

during the periodic motion. One can easily see from Fig. 2 that in the limit F — > the encircled area is equal to that of 
a triangle CAA'. With CB ~ Uq/F, AA' ~ 2Uq/F, we obtain the area A = (Uq/F) 2 and the Euclidean action is given 
by the expression S (0) = a(Uo/F) 2 (in Ref.i the result S (0) = (23/5) a(Uo/F) 2 has been obtained, in contradiction 
with the above analysis) . In the next section we shall generalize the zero temperature result to the case of arbitrary 
temperatures. 
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Fig. 2 Quasiclassical (x, j/)-trajectories corresponding to the tunneling of the vortex. The Euclidean action corresponding to the zero 
temperature trajectory (thick line) is equal to the encircled area. The thick dotted line marks a finite temperature bounce trajectory.. 



III. DECAY RATE 



A. Euclidean action 



As has been shown by VolovikEl (see also Refs. , the motion of a massless particle in a magnetic field subject 
to a potential U (x,y) is equivalent to the ID dynamics of a particle described by the Hamiltonian U(x,p/a). 
Consequently, our original 2D problem (|3|) reduces to the ID problem with the Hamiltonian given by the expression 



H(x,p) = U ^x 2 + (p/a) 2 ^j - Fx - 



(6) 



Let us shnw that in the limit F — > the problem is semiclassical. The semiclassical wave function can be written in 
the formEd 



T / \ C 

W (x) = —= exp 

*J x 



p dx \ — C' exp 



dx 



ip 1 d 2 H /dp 2 dp 
~h ~ 2 dH/dp 



(7) 



with x' an integration constant. The semiclassical approximation is applicable if 

d 2 H /dp 2 dp 



\p\ > ft 



dH / dp dx 



(8) 



Using Eq. (Iff) we can write for d 2 H/d p H 



^2 rr / a 2 



d 2 H/d P 2 



dH/dp 



1 

¥\ 



(p/a) 2 d 2 U /du 2 



(p/c 



(p/a) 



2 dllo/du 



(9) 



(2\ ^^ 2 2 
x 2 + (p/a) J .If x 2 + (p/a) <C a 2 , we can use Uq = nu 2 /2, i.e., \d 2 H/d p H\ = \/\p\ in the above 



limit. If 



(P/c 



> a 2 , 



%U /d p U ~ (x 2 + (p/a) 2 ) 



-1/2 



(10) 



a , we obtain 



and we obtain |<9 2 ir/<9piJ | < A/[p|, with A a constant of order one. Finally, at x 2 + (p/a)' 

again \d 2 H/d p H\ ~ A' ~ 1. Consequently, a sufficient criterion for the applicability of the semiclassical 

approximation takes the form 



h dp 
p 2 dx 



< 1, 



(11) 



which is the same criterion as for a "usual" Hamiltonian of the form H (x,p) = p 2 /2m + U (x) . 

Using the standard technique of binding senaiplassical wave-functions we obtain the following expression for the 
imaginary parts of the metastable energy levelsEl, 



r„ = 



u(E n ) 
An 



exp I ~ J \p\ dx J = 



u(E n ) 
An 



ex P I — r" 



(12) 
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with E n the energy levels at zero driving force, lu(E„) are the oscillation frequencies, and c„ and b n denote the turning 
points. The decay rate F can be found by averaging over the Boltzrnann distribution (a general discussion concerning 
finite temperature decay and the role of dissipation is found in Ref.tJ) 

r = (2/z)^r„ e - iS "/ T , (13) 

n 

with Z the partition function for the case F — 0. As the semiclassical approximation is applicable, we can substitute 
the sum in Eq. (|l3| ) by an integral and make use of the method of steepest descent. The extremal equation then 
takes the form dS/dE = —t(E), with t(E) the imaginary time oscillation period. Consequently, if one can calculate 
the function t(E) from the solution of the classical equation of motion, the function S(E) can be reconstructed via 

E 

simple integration, S(E) = — J T(E')dE' . Let us carry out this program for the present problem. 

The semiclassical trajectories can be found as the solution of the equation H (x,p) — E, with E the energy, i.e., 



p(x) = ±ay/f 2 (E + Fx)-x 2 , (14) 

with / = Uq 1 the inverse function of the potential shape Uq. There is a region c<a:<6(6>c>0), where 
the function p(x) is purely imaginary (the ^-coordinates c and b are associated with the turning points C and B 
in Fig. 2). The equation f 2 (E + Fx) - x 2 = has two solutions: At small x, f 2 (E + Fx) = (2/k) (E + Fx) , i.e., 
c = F/k + ^/F 2 /k 2 + 2E/k. Ash (Uq - E) /F , f 2 (E + Fx) ->■ oo, the equation f 2 (E + Fx) - x 2 = has another 
root b = (Uq — E)/F (the expressions for c and b are applicable for any energy not too close to Uq). The decay rate 



(2/h)f b c $\p\dx 



Note that everywhere inside the 



of a metastable state with an energy E is proportional to exp 

interval c < x < 6, except for the vicinity of the points c and b, \p\ = ax and the condition ([ll]) is fulfilled. At 
the points b and c, p = 0. These points play the role of turning points in the "usual" semiclassical approximationcj. 
Consequently, we have shown that in the limit F — > the semiclassical approximation is applicable. 
Let us calculate the Euclidean action. Using Eqs. (O) and (114) we can write for S 





S(E) = 2a J yjx 2 - f 2 (E + Fx) dx. (15) 



The oscillation time t(E) satisfies the equation t(E) = —dS/dE, i.e., 

r(E) = ^[ - Xdx (16) 
F J ^Jx 2 -f 2 (E + Fx) 

c 

In the limit F — > almost everywhere inside the interval [c, b] we have x 2 3> f 2 (E + Fx) and one can write for the 
period t(E) 

2a 

T (E) = — (b-c) + C(c) + C(b), (17) 

with C(c) and C(b) the contributions of the turning points c and b where the function x 2 — f 2 (E + Fx) vanishes. It 
can be shown that C(c) is relevant only if E is close to —F 2 /2k (see below) and the contribution of the point b is 
always negligible. Calculating C(c), substituting the result into Eq. (p^), and taking into account that for E not very 
close to Uq, b — c ~ (Uq — E)/F, we obtain 

r(E) ^^ {U Q-E) + ^\n f , (18) 



where x is an _E-independent cutoff parameter arising from the integration in the vicinity of the point c. 

Next, we need to solve the equation t(E) — h/T and find the energy E of the saddle-point trajectory at finite 
temperature T. The second term in Eq. (Qq) is relevant only if the solution of the equation t(E) — h/T is very close 
to —F 2 /2k, i.e., at low temperature where the equation 2a (Uq — E)/F 2 = h/T has no solution (this is the case when 
T < hF 2 /2aUo), and produces merely an exponentially small correction to the zero temperature result. This behavior 
is typical for a Hamiltonian problem where the finite temperature boundary conditions have a vanishingly small effect 
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on the bounce solution at small temperatures^. On the other hand, if T > hF 2 /2aUo, the solution of Eq. ([l8]) is given 

E 

by E = Uo — hF 2 /2aUo (up to exponentially small corrections). After a simple integration S(E) = — J r(E')dE' we 
obtain the Euclidean action SeucI = a(Uo — E) 2 / F 2 + E/T (the integration constant is obtained from the condition 
5(0) = alio 2 / F 2 ). In summary, for temperatures T < hF 2 /2aUo the Euclidean action is constant up to exponentially 
small corrections. At T\ — hF 2 /2aU , SeucI begins to decrease, SeucI (T) = hU /T — h 2 F 2 /4aT 2 . We then can write 
the following expression for the Euclidean action in the whole temperature interval (about the applicability of this 
result to the high temperature regime see below) 

T < hF 2 /2aU = Ti, 

(19) 

T > T t . 



Seu 



a 

frUo 
T 



h 2 F 2 
AaT 2 ' 



B. Crossover to Classical Behavior 



Next, let us calculate the crossover temperature T c from the thermal-assisted quantum regime to the purely thermal 
activation. Below we use the perturbative procedure which is applicable only for second-order transitions from 
quantum to classical behavior. For a first-order transition this approach breaks down. However, we will show that if 
the potential Uq(t) satisfies the required conditions (U(r) is monotonously increasing and U(r) ~Uq — B/r 2 , r — > oo) 
a second-order transition takes place. The crossover temperature T c is equal to Ti/tq, where To is the imaginary time 
oscillation period of the system in the vicinity of the time-independent thermal saddle-point solution. This solution 
is given by the equation x(t) = x max with x m ax the point where the function U(x,0) takes its maximal value. Near 



this point U(x,y) 
exactly, 



Uq — Bf {x 2 + y 2 ) — Fx. For this dependence the function t{E) (see Eq. (fL6|)) can be calculated 



t(E) 



2aB 



F 2 (b + d) 3/2 (b + c)(b-d) 1/2 



d)F 



+ (c - d) n 



tt (b-c)(b + c) (b-c)(b + c) 



2' (b-d)(b + d)']j {b-d){b + d) 
where b > c > d are the three roots of the equation 



Uo 



D 



Fx = E. 



7T (b-c)(b + C ) 

2 , V(b-d)(b + d) 



(20) 



(21) 



Again b and c are the turning points of the imaginary time trajectory, 

rir/2 



n(7r/2,n, k) 



is the complete elliptic integral of third order and 



(l — n sin 2 i 



1 — k 2 sir 



F(tt/2,A) 



7r/2 



1 - k 2 cos 2 



-1/2 



(22) 



(23) 



is the complete elliptic integral of second order. Here we consider only the physically relevant case E < E max = 
Uo - (3/2 2 / 3 ) B 1 / 3 F 2 / 3 . If E = E max , b = c = x max = (2B/F) 1/3 , and d = -l/2 2 / 3 (B / F) 1/3 . Substituting these 
values into Eq. ( po|) we obtain 



V3 fiF A / 3 
24/3^ aBi/a ■ 



(24) 



Inserting T > T c into Eq. (|19|) in the limit F — > we obtain Seuc\ (T > T c ) ~ Uq/T, such that we can use the result 
( |l9| ) at any temperature in this limit. The dependence Seu C i(T) is plotted in Fig. 3. 
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Let us show that the problem exhibits a second-order like transition at T — T c . The ID Hamiltonian system prnduces 
a smooth transition at T c if its imaginary time oscillation period t(E) is a monotonous function of energyEij E£l. A 
simple criterion to verify the monotonicity of the function t(E) is given by the derivative 8et evaluated at E max 
(see also RefES): For Be^e^* < the function is monotonous and we have a second-order like transition. If 
b ~ (Uq — E) 1 1 F ^> c, (i.e., for small F and E not too close to E max ), one can use Eq. ( |l7j ) to show that the function 
t(E) is monotonous in this interval. On the other hand, one can use Eq. (^p| ) as long as the roots b and c greatly exceed 

1/3 

the characteristic radius a of the pinning potential. In the limit F — > and for E close to E max , b, c ~ (B/F) 1 , and 
the condition 6, c S> a is well satisfied. 



>Eucl 



thermal assisted 
quantum 








2 nd order 



transition 



classical 



Fig. 3 Euclidean action as a function of temperature. At T < T\ = TiF 2 /2cxUq, SeucI is a constant up to exponentially small corrections. 
In the regime T > T\ the Euclidean action begins to decrease, see Eq. ([li]). The temperature T = T c (see Eq. (p4[)) marks the second-order 
like transition from quantum to classical behavior. In the vicinity of T c Eq. (^ij) gives an accurate description of the action. 



We reparametrize the energy in the form 



E = U - \B 1/3 F 2/3 



(25) 



where A = X(E) is a dimensionless parameter. Substituting this expression into Eq. ( pl| ) we obtain for large A (E 
away from E max ), c ~ X' 1 / 2 ^ / F) 1/3 , and b ~ X(B/F) 1/3 , i.e., 6 > c for A > 5 and Eq. © is valid in this region. 
Consequently, if we show that the oscillation period given by Eq. (OH) is a monotonous function of A for A < 5, a 

— 1 /3 

second-order transition from quantum to classical behavior takes place. The function (f = (rF /2a) (F/B) ) 



f(A) 



A 1/2 



(b + d)F 




( c -J)nf^,|- 5)(6 + 5) 



(&-£)(& + 2) 



2 (6-d)(6 + d) V (b-d)(b + d) 



(b-c)(b + c) 
(b-d)(b + d) 



(26) 



is plotted in Fig. 4 (solid line); here b > c > d are the roots of the equation 1/x 2 + x = X. We can see that the function 
f (A) is monotonously increasing. Surprisingly, even for A = A m i n ~ 3/2 2 / 3 , the slope of f is close to its asymptotic 
value f (A) ~ A + Const, A — > oo. Close to the point X m in = 3/2 2 / 3 , we find 



f(A) 



2V3 4V3 

^7T H 

V3 27 



7T A 



3 

22/3 



(27) 



and the derivative (?ET|£ max is indeed negative 
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dr/dE 



8\/37r aB 1 / 3 

F 4/3 



(28) 



After solving the equation t(E) = Ti/T in the vicinity of the point E = E max we obtain the expression for the action 

ti i . . ■"> -.1 ■' ,-.-■> ■■ \ .[■■—■• t-> _ ,2, 



(29) 



which improves on the result ( |l9| ) in the vicinity of T c 
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Fig. 4 The function f (solid line, see Eq. (|26[)), which is directly proportional to the imaginary time oscillation period t(E), as a 
function of A/A m ; n (see Eq. J25|)), A m ; n = 3/2 2 / 3 . Surprisingly, the slope of f near A m j n is very close to its asymptotic value at large A. 
The dotted line shows the function g(\) = + ^T"^ — ^7t) i see Eq. (p^). The dashed line illustrates the slope of the function 

/(A) = A + Const corresponding to the asymptotic expression for t(E), t(E) ~ 2a (U - E)/F 2 . Obviously, the slope of /(A) is very close 
to those of f(A) and g(A), indicating that one can use Eq. (IL8) in almost the whole energy range. 



C. Preexponential factor 

Finally, let us estimate the preexponential factor at high and at low temperatures. At T > T c the saddle-point 
solution is time independent and is given by the equation x = x max with x max the maximum of the function H(x, 0). 
In the vicinity of the point (x max ,0), 

H(x,p)~Uo 1 [X - Xmax) H o~ J V- (30) 

^max " ~^max 

The imaginary time Hamiltonian corresponding to (^0|) has the form of a harmonic oscillator with the mass a 2 x max /2B 
and the stiffness &^f-f max - Consequently at high temperatures one can use the result for the decay rate of a massive 
particle (see Refs.HO) 

sinh(^) / U \ 

hlI1 \^2 4 / 3 aB i / 3 T J 

This expression is .applicable at any temperature higher than T c except for a narrow temperature interval ~ ft 3 ^ 2 
around T c , see RefED. At T ^ T c , Tin/ a we obtain the simple result V = (hn/2-Ka) exp (—Uq/T). 

At zero temperature the decay involves only tunneling out of the ground state, i.e., T 

[cj(Eo)/2tt] exp [- (a/h) {U /Ff 



see RefS. Near the metastable minimum the vortex oscillates with the frequency 



= n/a, and we arrive at the result 
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27ra 



exp 



a /C7o 
F 



(32) 



However, one should point out that the preexponential factor in Eq. ( p2|) is only an order of magnitude estimate as 
the quasiclassical approximation is not properly applicable to the ground state. 



IV. APPLICABILITY 

Let us discuss the conditions for the applicability of the above results. First of all, we have to account for a nonzero 
vortex mass. In this case, the system is described by the Hamiltonian 



H = 



P' 2 



(Py 2 



ax) 



2m 



U, 



(v^ 



V 



- Fx, 



(33) 



which is equivalent to the Hamiltonian of a massive 2D charged particle in a magnetic field oriented perpendicular 
to the plane of motion. Let us discuss the conditions which should be satisfied in order to use the Hamiltonian (||) 
instead of that given by Eq. (33). It is useful to introduce the four new operators (see also Ref.El) £, y — ~l 2 P x /h, 
X = —l 2 P y /h, £ x = l 2 P y /h + x, Y = l 2 P x /Ti + y, where I = y/h/a is the magnetic lengthEl. We assume that the 



yj lb T O,, 1 - 1 il/HTj/, 

condition I <C a is satisfied (below we will verify this condition using parameters for high-X^-supcrconductors). These 

new operators satisfy the commutation relations = [X, Y] = il 2 , [<£ X ,Y] = = fe/i^] = = 0- 

Going over to the X, ^-representation the Hamiltonian ( |33| ) can be written in the form 



H 



2m d& 



2 >-2 



(34) 



In a classical language the problem we study involves two types of motion: A fast rotation of the particle with the 
cyclotron frequency lu c — a/m superimposed on the slow guiding center motion. The term 



h= 



n 2 d 2 



muJ cC 



2m dQ 



(35) 



describes the "fast" part of the Hamiltonian. The lowest eigenvalue of h is equal to fiuj c /2. If the characteristic 
variation of the potential on the length I is much smaller than fko c /2, we can average the Hamiltonian ([34]) over the 
groundstate eigenfunction x(£y) °f the operator h, see Refo. This process corresponds to the averaging over the fast 
rotation with the frequency u) c . After averaging we obtain 



H^(H) iast = / C/(k)exp 



k 2 l 2 



+ i (k x X + kyY^j 



d 2 k. 



(36) 



Note that in the X -representation Y = 
On the other hand, if I is small, [X, Y] - 
Hamiltonian 



~il 2 d/dX. If I -C a, we can neglect the exponent exp (— k 2 / 2 /4) in Eq. (|3f 
and we can easily perform the integration in Eq. ([56]) to arrive at the new 



H eS = U \^X 2 + Y 2 



FX 



(37) 



describing the guiding center motion. As the operators X and Y obey the commutation relation [X, Y] = il 2 , we can 
identify X with the coordinate and Y with the momentum. This Hamiltonian is identical to that in Eq. (^|). 

In order to carry out the above procedure the characteristic variation of the potential U(x, y) on the scale I has to be 
much smaller than Tko c , i. e., kI 2 -C huj c , and with k ~ Uo/a 2 we obtain the condition a 2 a 2 / mU^ 3> 1. Consequently, 
we arrive at the following two conditions for the applicability of the approach we have used, 



I <C a and 



mUr 



» 1. 



(38) 



Let us estimate the parameters in Eq 
2 



(pq) appropriate for high-T c -superconductors: a — irhns, where n = 



10 21 cm 3 is the electron density and s — 15 A is the interlayer spacing in Bi-2:2:1:2. For the magnetic length we 







then obtain I = \fhj~a. ~ 4 A <C a ~ 30 A. For the second condition in Eq. (38) we use the expression for the vortex 
mass predicted for the superclean limitEl m = m e (ei?/A) 2 ~ 100m e and the depth of the potential as given by the 
equation Uq = ($o/47rA a f,) s, with $o the magnetic flux quantum and X a b = 2000 A the a6-magnetic penetration 
length. Finally, we obtain a 2 a 2 /mUo ~ 10 ^ 1, i.e., weJind that both conditions in Eq. ( |38| ) are well satisfied. 

The above results should be compared to those in Ref .□ where the condition ha/mUo 3>i 1 for the applicability of the 
results has been used. Estimating this parameter we obtain ha/mU&.~ 0.1 <C 1. In Ref.EI an estimate Ha/mUo — 200 
has been obtained using the vortex mass from the dirta limitE^TO, which is much smaller than the mass in the 



superclean limit. We see that the LLL-approach of Ref.cl breaks down in the superclean limit if one does not take 
explicitly the condition I <C a into account. 

Finally, let us show that the second condition in Eq. ([38]) can be easily obtained from a simple analysis of the 
semiclassical bounce trajectories: If the Euclidean action is much larger than unity, we can use the instanton method 
for the calculation of the imaginary part of the partition function which determines the decay rate T cx ImZ /Z. In 
order to determine ImZ within exponential accuracy we have to solve the classical equation of motion describing the 
bounce solution. If it turns out that the correction of the bounce trajectory due to the mass term is small, we arrive at 
an effective Hall tunneling problem. It is possible to neglect the mass if everywhere along the trajectory m \v\ -C a\v\, 
where v is the characteristic acceleration of the vortex during the imaginary time motion and v its characteristic 
velocity. If a vortex is moving along the sides CA,AA', or A'C of the triangle CAA' (see Fig. 2), we can make the 
estimate |«| ~ v 2 /(Uo/F), whereas for v the estimate v ~ Uo/aa holds, i.e., the following condition should be satisfied 

F a 2 a 2 
b c mU 

with F c ~ Uo/a the critical force. On the other hand, in the vicinity of the point C, we can write v ~ Uo/aa and 



u ~ v 



2 /a, i.e., the condition a 2 a 2 /mUo ^> 1 has to be satisfied, which is identical to the second condition in Eq. (|3 



Comparing this result with Eq. pH) we see that the second condition in Eq. (38) is stronger since F < F c 



V. SUMMARY 



In section IV we have described how to reduce the dynamics of a massive particle in a magnetic field to the guiding 
center motion along equipotential lines of a smooth potential. Carrying out this procedure for vortices is not entirely 
unproblematic. In fact, we have seen, that the averaging process over the fast component introduces a frequency 
lu c = a/m ~ wo which is at the limit of the applicability of the low frequency vortex dynamics as described by 
Eq. ([!]). This problem has been ignored in previous works as the inertial term has been dropped on the classical level 
of Eq. As shown above, both methods produce the same condition a 2 a 2 /mUo 1 for the irrelevance of the mass 
term, however, starting from the description ( j33|) one obtains that the frequency lo c — a/m naturally shows up in the 
quantum description. 

Briefly summarizing, we have considered the decay of the metastable state of a massless particle in a magnetic 
field trapped in a cylindrical attractive potential and subject to a small external force. The above problem is an 
appropriate model for the depinning of pancake vortices in superclean Josephson-coupled superconductors for the 
case of small (J -C j c ) transport currents. The Euclidean action as a function of temperature is given by Eq. (|lS|); 
For T < T\ the action is constant up to exponentially small corrections and decreases smoothly above T\. At T = T c , 
the crossover to the calssical regime takes place. Close to and above T c the decay rate can be accurately described 
by Eq. (^9|) which reduces to the exact classical result above T c . For the potential satisfying the conditions i) U{r) is 
a monotonous function of r and ii) U(r) = Uq — B/r 2 , r —> oo, the transition from quantum to classical behavior is 



second-order like and the crossover temperature T c is proportional to F 4 / 3 , see Eq. (24) 
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